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Abstract. A pointed observation of the galactic high-mass X-ray binary 4U 0114+65 was carried out with 
BeppoSAX to compare the X-ray spectral and timing characteristics observed by this satellite over the broadest 
range of energies thus far (1.5-100 keV) with the information previously obtained with other spacecraft. The light 
curve of 4U 0114+65 shows a large flare at the beginning of the BeppoSAX pointing and no significant hardness 
evolution either during the flare or in the low state occurring after the flare itself. The modulation at ~2.7 hours, 
attributed to the accreting neutron star (NS) spin periodicity, is not significantly detected in our data, although 
fluctuations with timescales of -~3 hours can be seen in the 2-10 keV light curve. Shorter modulations down to 
timescales of minutes are also found and interpreted as due to accretion of matter onto the NS. The flaring and 
the low state spectra of 4U 0114+65 can be equally well fitted either with a power law modulated by a high-energy 
exponential cutoff or with a Comptonization model. During the low state the presence, although tentative, of a 
thermal component (with kT ~ 0.3 keV) at low energies, possibly produced by an ionized plasma cloud around 
the NS, cannot be excluded. Contrary to previous claims, a cyclotron resonant feature in absorption at ~22 keV 
was not detected in the BeppoSAX spectroscopic data, whereas evidence for a Fe emission line around 6.4 keV is 
found only during the low state emission. Using all of the above information, a scenario for the system in which 
the NS is embedded in, and accreting from, a low angular momentum gas cloud is envisaged. 
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1. Introduction 

Persistently emitting High-Mass X-ray Binaries (HMXBs) 
are gravitationally bound systems composed of a early- 
type star losing matter, generally via stellar wind, in 
favour of a compact object, a neutron star (NS) or possi- 
bly a black hole. This accretion phenomenon induces high- 
energy emission from the accreting flow and, in the case 
of a NS as the accretor, from the surface of the compact 
object. This kind of HMXB has X-ray luminosities of the 
order of at least 10 35 erg s _1 in the 2-10 keV range. Also, 
if the accreting object is a magnetized NS, the accretion 
flow is channeled onto the magnetic polar caps and the 
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X-ray emission can appear as pulsed, modulated by the 
NS spin period (see White et al. 1995 for a review). 

The hard X-ray source 4U 0114+65 belongs to this 
group of X-ray binaries. It was discovered by the SAS-3 
satellite (Dower et al. 1978) and extensively observed in 
the X-ray band by several spacecraft. 

The lightcurve of the object appears quite variable in 
X-rays, showing flares lasting a few hours, occurring on 
daily timescales and that are 15-20 times brighter than 
the persistent low-level emission (Yamauchi et al. 1990; 
Apparao et al. 1991). Shorter-term flickering on timescales 
of minutes was also observed (Koenigsberger et al. 1983). 
Finley et al. (1992), by analyzing archival EXO SAT and 
ROSAT data, showed the existence of a very stable 2.78- 
hour periodicity in the X-ray light curve of 4U 0114+65; 
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this was more recently confirmed by Corbet et al. (1999) 
and Hall et al. (2000) by using RXTE ASM and PCA 
data, respectively. This was interpreted by these authors 
in terms of the the spin period of a slowly rotating NS. 
Farrell et al. (2005) moreover reported the presence of a 
superorbital periodicity of 30.7 days in this system, again 
by the analysis of RXTE/ ASM data. This may suggest 
the presence of a precessing accretion disk around the NS. 

X-ray spectroscopy (Yamauchi et al. 1990; Apparao et 
al. 1991; Hall et al. 2000) shows that the source spectrum 
is satisfactorily fitted with the 'classical' phenomenological 
model for accreting NSs in HMXBs (White et al. 1983), 
i.e. a power law modified at high energies by an expo- 
nential cutoff. In the case of 4U 0114+65, typical spectral 
parameter values are found to be T ~ 1, E cut ~ 8 keV and 
Efoid ~ 20 keV for the photon index, the cutoff and folding 
energies. A neutral iron emission line at 6.4 keV was also 
sometimes detected, preferentially during the persistent 
low state emission. 

Bonning & Falanga (2005) observed this source with 
INTEGRAL and again found the 2.7-hour spin mod- 
ulation up to an energy of 80 keV. The 5-100 keV 
INTEGRAL spectrum was consistent with the findings, 
illustrated above, of previous high-energy observations. 
They also measured a NS spin-up trend P = — 8.9xl0 -7 
s s" 1 over ~8 years. Moreover, these authors tentatively 
identified a Cyclotron Resonant Feature (CRF) at ~22 
keV in the X-ray spectrum of the source, implying a mag- 
netic field of 2.5xl0 12 G for the accreting NS. 

This evidence can be interpreted as due to accretion 
onto a NS from a wind coming from an early-type star. 
Indeed, 4U 0114+65 was associated with its optical coun- 
terpart, the 11 th magnitude star V662 Cas (also known as 
LSI +65°010) by Margon & Bradt (1977). Subsequently, 
Crampton et al. (1985) spectroscopically determined the 
orbital period of the system as P or b = 11^59; this modula- 
tion was also found later in the RXTE/ ASM X-ray data 
by Corbet et al. (1999). Reig et al. (1996) performed a 
thorough study of the optical properties of the companion 
star and found it is a blue supergiant of spectral type Blla 
located at 7.0+3.6 kpc from earth, though Koenigsberger 
et al. (2003) note that it could be located closer. This 
star is in a basically circular orbit (Crampton et al. 1985; 
Koenigsberger et al. 2003) and may induce X-ray eclipses 
(Hall et al. 2000). 

The above optical information conclusively placed 4U 
0114+65 among wind-accreting supergiant HMXBs in 
terms of spectral type of the companion and of X-ray 
luminosity budget, the latter being between 10 35 and 10 36 
erg s _1 in the 2-10 keV range; indeed, previous distance 
estimates (d ~ 2 kpc; e.g., Aab et al. 1983) led to a 2- 
10 keV luminosity of about ten times lower, erroneously 
implying that this source was a low-luminosity HMXB. 

No radio emission was detected from this source in the 
5 GHz band down to a level of 1.2 mJy (Nelson & Spencer 
1988). 

We observed 4U 0114+65 with BeppoSAX (Boella et 
al. 1997a) within the framework of our program of pointed 



observations of HMXBs (Dal Fiume et al. 2000; Orlandini 
et al. 1998,1999,2000,2004; Masetti et al. 2004) in order to 
describe its X-ray temporal and spectral behaviour over 
the wide range of energies (0.1-300 keV) covered by this 
satellite, and to compare its X-ray characteristics with 
those of other objects belonging to this class. 

Here we present the results of a pointed observation 
made with BeppoSAX on this system in January 1999. 
Thanks to the broadband spectroscopic capabilities of this 
satellite and to its high sensitivity, for the first time we 
were able to simultaneously explore the 1.5-100 keV emis- 
sion from this source. The paper is organized as follows: 
Sect. 2 describes the observations and the data analysis, 
in Sect. 3 the results showing the X-ray spectral and tim- 
ing behaviour of 4U 0114+65 will be reported; in Sect. 4 
a discussion is given. 

2. The BeppoSAX observation 

4U 0114+65 was observed between January 26 and 
27, 1999, for ^63 ks, by the four coaligned Narrow- 
Field Instruments (NFIs) carried by BeppoSAX: the Low 
Energy Concentrator Spectrometer (LECS, 0.1-10 keV; 
Parmar et al. 1997), two Medium Energy Concentrator 
Spectrometers (MECS, 1.5-10 keV; Boella et al. 1997), 
the High-Pressure Gas Scintillation Proportional Counter 
(HPGSPC, 6-60 keV; Manzo et al. 1997) and the 
Phoswich Detection System (PDS, 15-300 keV; Frontera 
et al. 1997). The total duration of this BeppoSAX point- 
ing along with the on-source exposure times for each NFI 
used are reported in Table 1. 

The source was well detected by all four NFIs. Good 
NFI data were selected from intervals outside the South 
Atlantic Geomagnetic Anomaly when the elevation angle 
above the earth limb was >5°, when the instrument func- 
tioning was nominal and, for LECS events, during space- 
craft night time. The SAXDAS 2.0.0 data analysis package 
(Lammcrs 1997) was used for the extraction and the pro- 
cessing of LECS and MECS data. The PDS data reduction 
was performed using XAS version 2.1 (Chiappetti & Dal 
Fiume 1997). LECS and MECS data were both reduced 
using an extraction radius of 4' centered on the source po- 
sition; before extraction, data from the two MECS units 
were merged. Background subtraction for the two imaging 
instruments was performed using standard library files. 
The background for the PDS data was evaluated from 
the fields observed during off-source pointing intervals, 
whereas that of the HPGSPC was computed through an 
Earth-occultation technique (Manzo et al. 1997). 

Because 4U 0114+65 is located near the Galactic plane 
and is not a particularly bright X-ray binary, we checked 
for possible effects induced by Galactic diffuse emission 
in the PDS data background evaluation. The PDS off- 
source fields for background evaluation were indeed at dif- 
ferent Galactic latitudes {b = -0?9 and b = +6?0) with 
respect to the source (b = +2?6), so a possible gradient in 
the Galactic diffuse emission, if present, can be measured 
and eventually removed from the PDS data. However, the 
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Fig. 1. Background-subtracted X-ray light curves of 4U 0114+65 in the 2-5 keV (upper panel), 5-10 keV (middle 
panel) bands, along with their ratio (lower panel) as observed during the BeppoSAX pointing. All curves are rebinncd 
at 300 s. A large flare during the first 3 hours of the observation is apparent. Subsequent flickering activity down to 
timescales of minutes can be noticed. Time is measured from the start of the BeppoSAX observation, on January 26, 
1999, at 21:38:42 UT (see also Table 1), and it is expressed in hours starting from 00:00 UT of January 26, 1999. 
Vertical error bars show la confidence level uncertainties for each bin. 



count rate difference between the two fields is — 0.06±0.07 
counts s , thus consistent with zero; moreover, this dif- 
ference impacts on the background estimate by less than 
0.3%. Therefore we considered this effect negligible. 

3. Results 

3.1. Light curves 

We extracted the 2-5 and 5-10 keV light curves of 

4U 0114+65 from the BeppoSAX/MECS data. These 



background-subtracted MECS light curves rebinned at 
300 s, along with their ratio, are shown in Fig. 1. A 
large flare, similar to those reported by Apparao et al. 
(1991), was seen from 4U 0114+65 at the beginning of 
the BeppoSAX observation (Fig. 1). Its peak intensity is 
~20 times the persistent low-level emission observed after 
the event. No significant hardness evolution was detected 
during the episode, although the (5-10 keV)/(2-5 keV) 
hardness ratio seems to increase starting ^5 hours after 
the beginning of the observation, just after the end of the 
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2-10 keV intensity (Cts/s) 

Fig. 2. Hardness-intensity diagram for the 2-10 keV emis- 
sion of 4U 0114+65 during the BeppoSAX pointing. 
Although an hint of an inverse correlation between the 
5-10 keV / 2-5 keV hardness ratio and the 2-10 keV in- 
tensity is suggested by the plot, no statistically significant 
trend is present. Error bars show la confidence level un- 
certainties for each data point. 

flare (Fig. 1, lower panel). The light curves also show sub- 
stantial random variability in the form of internal fluctu- 
ations lasting from a few hours to minutes (the so-called 
'X-ray flickering'). 

In order to see whether this erratic variability implied 
spectral changes depending on the source intensity, we 
computed a hardness-intensity ratio between the 5-10 keV 
and the 2-5 keV count rates and plotted it against the to- 
tal 2-10 keV count rate. The results are plotted in Fig. 2: 
although the figure hints to an inverse correlation between 
the two quantities, there is no significant dependence of 
the hardness ratio on the total intensity of 4U 0114+65 
in the 2-10 keV range, as opposed to the inverse relation 
between X-ray hardness and intensity found by Apparao 
et al. (1991) from EXOSAT data.. 

Using the optical orbital ephemeris reported by 
Crampton et al. (1985) it is found that this large flare oc- 
curred just after periastron passage: in detail, we find that 
the flare occurred at phase <p ~ 0.05 assuming the most 
likely case of circular orbits for the system (Koenigsberger 
et al. 2003). 

3.2. Timing properties 

Timing analysis on the 2-10 keV data (where the S/N ra- 
tio was highest) was performed with the FTOOLS V5.1 1 
(Blackburn 1995) tasks powspec and efold, after hav- 

1 available at: 
http : / /heasarc . gsfc .nasa. gov/f tools/ 
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Fig. 3. 2-10 keV PSD of 4U 0114+65 obtained from the 
events registered by the two MECS units. The PSD is 
normalized according to the prescription by Leahy et al. 
(1983). No coherent oscillation or periodicity is found, 
whereas a shot- noise l//-type trend is apparent, flatten- 
ing at frequencies below / w 10~ 4 Hz. 

ing converted the event arrival times to the solar system 
barycentric frame. 

The results do not reveal the presence of any kind of 
periodicity or quasiperiodic oscillations. We were not able 
to detect the 2.7 hour NS spin periodicity in our analysis. 
The Power Spectral Density (PSD; Fig. 3) obtained with 
this analysis is characterized by red noise and shows no sig- 
nificant deviations from the l//-type distribution (where 
/ is the time frequency) typical of the 'shot-noise' be- 
haviour. This trend however flattens for frequencies lower 
than / w 10~ 4 Hz, which is of the order of magnitude of 
the 4U 0114+65 NS spin frequency. Indeed, one can see 
from Fig. 1 that, at least in the first half of the BeppoSAX 
observation, the presence of fluctuations on a timescale of 
^3 hours is suggested. We stress that the non-detection 
of this modulation may be due to the relatively short du- 
ration of the observation (~3 pulse cycles if one considers 
the on-source time), to the sparse light curve sampling 
and to the presence of the large flare at the beginning of 
the pointing. 

3.3. Spectra 

To perform the spectral analysis, the NFI pulse-height 
spectra were rebinned to oversample by a factor of 3 the 
full width at half maximum (FWHM) of the energy reso- 
lution and to have a minimum of 20 counts per bin, such 
that x 2 statistics could reliably be used. Data were then 
selected, for each NFI, in the energy intervals in which 
sufficient counts were detected from the source and for 
which the instrument response function was well deter- 
mined. This led us to consider the spectral interval 1.5-4 
keV for the LECS, 1.6-10 keV for the MECS, 7-30 keV 
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Fig. 4. Left panel: 1.5-100 keV X-ray spectrum of 4U 0114+65 obtained with the BeppoSAX NFIs during the flare 
reported in Fig. 1 and fitted with a photoelectrically absorbed PL modified with a high energy exponential cutoff 
(HighECtjtPL model of Table 2). The best-fit model is shown as a continuous line. (Right panel): 1.5-100 X-ray 
spectrum of 4U 0114+65 obtained with the BeppoSAX NFIs during the low-level state following the flare reported in 
Fig. 1. The best-fit model is the same as in the left panel. Note the hint of a soft X-ray excess below 3 keV. For both 
panels, the residuals between data and model are shown in units of a, and the differences in counts s _1 keV -1 among 
the various instruments arc due to different effective detector areas. 



for the HPGSPC and 15-100 keV for the PDS. We then 
used the package XSPEC (Dorman & Arnaud 2001) vll.3.1 
to fit the resulting broad band energy spectrum. 

We included in all fits described here an inter- 
stellar photoelectric absorption column, modeled using 
the Wisconsin cross sections as implemented in XSPEC 
(Morrison & McCammon 1983) and with solar abundances 
as given by Anders & Grevesse (1989). 

When performing the spectral fits, normalization fac- 
tors were applied to LECS, HPGSPC and PDS spectra 
following the cross-calibration tests between these instru- 
ments and the MECS (Fiore et al. 1999). These factors 
were constrained to be within the allowed ranges dur- 
ing the spectral fitting, with the possible exception of the 
HPGSPC one, with results slightly (-10%) higher than 
the expected range. This shift is however justified by the 
fact that the MECS and PDS data, in the ranges in which 
they overlap those of the HPGSPC, follow the same bin- 
to-bin variation pattern. We can moreover exclude the 
presence of a contaminating field source because (i) no 
catalogued sources with 2-10 keV flux >10~ 13 erg cm~ 2 
s _1 , i.e. 10 3 times fainter than 4U 0114+65, are present in 
the HPGSPC field of view (l°xl° FWHM), (ii) the MECS 
and PDS instruments do not need such a higher shift for 
their intercalibration factor and, most importantly, (iii) 
no source other than 4U 0114+65 is detected in the field 
imaged by the MECS (which has a size comparable to that 
of the HPGSPC). 

To explore possible spectral differences between the 
flaring and the low-level emission from 4U 0114+65 which 
could be overlooked by the simple 2-10 keV hardness- 
intensity ratio analysis of Fig. 2, we divided the observa- 
tion into two parts using the data before and after 01:00 
UT of January 27, 1999 (i.e., the 25 th hour in the plot in 
Fig. 1). This allowed us to construct broad (1.5-100 keV) 



'flare' and 'low-level' post-flare spectra for this source (re- 
ported in Fig. 3), therefore making the analysis of possible 
spectral changes described in Sect. 3.1 more complete. For 
these two time-resolved spectra Table 2 reports the mod- 
els along with the corresponding parameters that best fit 
the data. 

In the flaring spectrum, the 'classical' phenomenolog- 
ical model of a power law (PL) plus a high-energy ex- 
ponential cutoff (HighECtjtPL; White et al. 1983) pro- 
vides an acceptable fit to the data, with best-fit parame- 
ters which are similar to, but more accurately determined 
than those found e.g. by Bonning & Falanga (2005) on 
an INTEGRAL X-ray spectrum with similar broadband 
coverage. An equivalently good fit (see Table 2) was also 
achieved using a PL with a simple cutoff (CutoffPL), in 
the form N(E) oc • e ~ E l E ™K 

Turning to a more physical description, as suggested 
by recent observational (e.g., Masetti et al. 2004; Torrejon 
et al. 2004) and theoretical (Becker & Wolff 2005) results, 
we tried a Comptonization description of the source flare 
spectrum via the model by Titarchuk (1994; COMpTT in 
xspec) and found that it also provided an acceptable de- 
scription, with best-fit parameters reported in Table 2. We 
found no statistically significant difference in the fit when 
using different geometries for the Comptonizing cloud. We 
therefore decided to assume the (default) spherical geom- 
etry. Also, using the simpler, non-rclativistic modelization 
by Sunyaev & Titarchuk (1980; COMpST in xspec) for the 
Comptonization, we obtained acceptable but nevertheless 
worse fits (see Table 2); in particular, this description sys- 
tematically underestimates the data of the flare spectrum 
above 40 keV. 

All of these models also describe the low-level spectrum 
of 4U 0114+65 above 3 keV acceptably well. However, 
as one can see in Fig. 4, below this energy an apparent 
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Table 2. Best-fit spectral parameters for the flaring and 
low-level emission of 4U 0114+65. In the leftmost column 
of the table, ATjj is the neutral hydrogen column density, 
T is the power law photon index, E cut is the cutoff energy, 
Efoid is the exponential folding energy, T is the temper- 
ature of the Comptonized seed photons, T c - is the tem- 
perature of the Comptonizing electron plasma, r is the 
optical depth of the Comptonizing plasma, and the Ks 
are the normalization constants of each model. Quoted 
errors are at the 90% confidence level for a single param- 
eter. Luminosities, corrected for interstellar Galactic ab- 
sorption, are computed assuming a distance d = 7.0 kpc 
and using the HighECutPL model as reference; no sig- 
nificant differences arc found if other models are used. 



Model + parameters 


Flare 


Low-level 


CutoffPL: 

r 

Scut (keV) 
K PL (xlO" 2 ) 
XVdof 


s Q+°- 8 
0.86+0.15 

9 9 +0.6 
^• z -0.4 

201.7/199 


17.61?;? 
q9+0' ±9 
18l| 

o-7±g:i 

273.2/194 


HighEcutPL: 

N H a 

r 

Scut (keV) 
S fo i d (keV) 
K PL (xlO" 2 ) 
X 2 /dof 


q 7+0.7 

y -'-0.9 

■t qq + 0.09 

12+ 2 
2ltt 
3.4l°oi 
199.0/198 


15.41?;? 

0.9+0.2 
6.0l° ;? 
17l 5 3 
0.43l°;?i 
262.7/193 


CompTT: 
N H a 

To (keV) 
T e _ (keV) 
r 

-KCompTT (Xl0~ 3 ) 

XVdof 


cc c+0.9 
Q.O_q 7 

-| Q1+0.12 

1 - J1 -0.14 

4.il° :S 

212.7/198 


H-7lii 
1.51+0.15 

q 7+2-4 

9 - ' -1.6 

3.5+0.6 
1.4+0.3 
261.9/192 


CompST: 
To (keV) 

T 

^CompST (Xl0~ 2 ) 

XVdof 


10.7+0.7 

°- i -0.5 

10.4+0.9 

c 9+0.8 

229.9/199 


20.9+1.5 

'■°-1.0 

8.8+?? 
1 q+0-4 

290.0/194 


T b 

^2-10 keV 

r b 

MO-100 kcV 


13.5 
24.1 


3.7 
7.9 


a In units of 10 22 cm" 2 
b In units of 10 35 erg s" 1 



soft X-ray excess is found despite the high value of iV H , 
with the overall 1.5-100 keV fit providing a x 2 /dof ~ 1.3- 
1.4 (where 'dof stands for 'degrees-of- freedom') indepen- 
dently of the chosen model among those listed in Table 2. 
We tried to fit this excess with a thermal component, and 
found that the fit is marginally improved (to a reduced 
X 2 ~ 1.2, which corresponds to a 90% confidence level ac- 



cording to the F-test; e.g., Press et al. 1992) using either 
a blackbody (BB) with kT t h ~ 0.3 keV and radius R t h ~ 
200 km, or with a disk-BB model (DBB; Mitsuda et al. 
1984) with similar temperature and radius, or with a hot 
diffuse thermal plasma (mekal model in xspec; Mewe 
et al. 1985), again with a similar temperature as of the 
BB model. In the case in which this thermal component is 
added to the Comptonization, an acceptable fit (x 2 /dof = 
237.6/191) is also obtained if one forces the thermal emis- 
sion to be the source of the Comptonization seed photons 
by imposing kT t h = kT . The value of this parameter is 
again ^0.3 keV, whereas those of the Comptonization do 
not differ from those reported in Table 2. 

In Table 2 the unabsorbed luminosities of the source 
in the two brightness states and in the 2-10 and 10-100 
keV ranges are also reported, assuming a distance to 4U 
0114+65 of 7.0 kpc (Reig et al. 1996). To compute them, 
the best-fit PL plus high energy exponential cutoff (i.e., 
HighECutPL) modeling described above as the reference 
spectral modelization of the source was used. No substan- 
tial differences in these luminosity values are found if one 
uses the PL plus simple cutoff or the Comptonization best- 
fit models. 

Finally, we searched for the presence of an iron emis- 
sion line around 6.4 keV (modeled as a Gaussian) and 
of a CRF at ~22 keV (using the CYCLABS multiplica- 
tive model in xspec ; Mihara et al. 1990; Makishima et 
al. 1990). From our best-fit descriptions there appears to 
be no need to include either of the two spectral features 
mentioned above in the 4U 0114+65 flaring state. Instead, 
concerning the low-level state, while we again did not find 
any evidence for a CRF, a broad (FWHM ~ 1 keV) Fe 
emission is required with a probability of improvement by 
chance of 5.8 xl0~ 5 . 

For each of the two parts of the BeppoSAX observa- 
tion we thus computed the value (when applicable) or the 
90% confidence level upper limits to the equivalent width 
(EW) of a Fe emission at 6.4 keV (assuming different line 
FWHMs) and to the depth of a CRF at 22 keV (assuming 
a line width of 10 keV, as found by Bonning & Falanga 
(2005) and which is typical of magnetic NSs in HMXBs; 
e.g. Orlandini & Dal Fiume 2001). The results, reported 
in Table 3, are independent of the chosen best-fit spectral 
description. 

However, from a visual inspection of the right panel of 
Fig. 4, there is a hint of a feature around 20 keV in the low- 
level spectrum of 4U 0114+65. In order to check whether 
the feature is real or due to a non-perfect modeling of the 
continuum we performed a Crab ratio of the PDS spec- 
tra for both brightness states. This technique was used 
to successfully pinpoint CRFs in numerous X-ray binary 
pulsars (see, e.g., Orlandini et al. 1998), and consists of 
obtaining the ratio between the count rate spectrum of 
the source and that of the featureless power law spectrum 
of the Crab. This ratio has the advantage of minimizing 
the effects due to detector response and calibration un- 
certainties. Any real feature present in the source spec- 
trum can then be enhanced by multiplying the ratio by 
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Table 3. Values and 90% confidence level upper limits to 
the EW (in eV) of the 6.4 keV Fe emission line, assuming 
three different fixed line FWHMs, and to the CRF depth 
assuming an energy Eqrf = 22 keV and a width of 10 
keV. 



Observation 


Fe 


line EW (eV) 


CRF depth 




0.1 keV 


0.5 keV 1 keV 




Flare 


<61 


<93 <141 


<0.3 


Low state 


80±50 


250±go° 580+200 


<0.3 



the functional form of the Crab spectrum (a power law 
with photon index T = 2.1) and by dividing it by the 
functional form of the continuum adopted to fit the broad- 
band source spectrum, obtaining the so-called Normalized 
Crab Ratio (NCR; see Orlandini et al. 1998 for details). 
The NCR performed on both the flare and low-level 4U 
0114+65 spectra docs not show any incorrect modeling of 
the continuum and further excludes the presence of any 
absorption feature at ~22 keV. 

Also, if we force our low-level spectral data to be de- 
scribed by a spectrum like that of the best fit by Bonning 
& Falanga (2005), we obtain a larger reduced x 2 (~l-8), 
a larger HPGSPC intercalibration constant, and a wide 
fictituous emission bump between 20 and 50 keV, i.e., a 
trend in the fit residuals which is opposite with respect to 
that found by those authors. Because of this, we regard 
the description of our BeppoSAX low-level data using the 
best-fit results by Bonning & Falanga (2005) as statisti- 
cally not satisfactory. 

4. Discussion 

The BeppoSAX observation of 4U 0114+65 presented in 
this paper allowed us to explore the X-ray behaviour of 
this source in the widest band thus far achieved (1.5-100 
keV). Analysis of the X-ray light curves of the object 
shows that it underwent a large flare, with total dura- 
tion ~4 hours, at the beginning of the pointing. This was 
followed by a progressive return of the source to a low- 
level flux at an intensity about 20 times fainter than the 
peak of the flare. One can indeed see from Table 2 that 
the average flux level during the flare is more than ^3 
times that of the average post-flare emission, going from 
3.8xl0 36 erg s" 1 down to 1.2xl0 36 erg s" 1 in the 2-100 
keV band. Using the orbital ephemeris determined opti- 
cally by Crampton et al. (1985) we found that this flare 
occurred just after periastron passage, at phase </> ~ 0.05. 

The modulation at ^2.7 hours, attributed to the ac- 
creting NS spin periodicity, is not significantly detected 
in our data, although fluctuations with timescales of ^3 
hours could be seen in the 2-10 keV light curve, and the 
PSD of the source in this X-ray range indicates a break 
around the corresponding temporal frequency (ftT0~ 4 



Hz). This may be due to the relatively short duration of 
the observation (~3 pulse cycles if one considers the on- 
source time), to the sparse light curve sampling and to the 
presence of the large flare at the beginning of the pointing. 
X-ray flickering down to timescales of minutes was instead 
found, with its PSD being proportional to 1//, where / 
is the temporal frequency. These short-term variations in 
the X-ray flux may be due to the instability of the ac- 
cretion process, or to the inhomogeneities in the accreting 
stellar wind (or both; see, e.g., Kaper et al. 1993). Thus, 
both of these long-term (flares) and short-term (flicker- 
ing) 'shot-noise' variabilities point to an explanation for 
this X-ray activity as due to random inhomogeneities in 
the accretion flow onto a compact object (e.g. van der Klis 
1995). 

The possible presence of spectral evolution was inves- 
tigated through a 2-10 keV hardness-intensity diagram 
and by accumulating the 1.5-100 keV X-ray spectra of 
the flare and of the low-level post-flare emission. The 
hardness-intensity diagram does not decisively show ap- 
preciable color evolution in the 2-5 keV and 5-10 keV 
bands flux; on the other hand, the comparison of the 
flare and low-level spectra of 4U 0114+65 as observed by 
BeppoSAX shows that (leaving aside the obvious varia- 
tion in the model normalization due to the change in the 
flux level) no substantial modifications in the spectral pa- 
rameters are apparent, with the notable exception of the 
neutral hydrogen column density Ah- 

Indeed, its value nearly doubles at lower emissions in- 
dependently of the chosen best-fit spectral modclization. 
This trend was already suggested by the EXOSAT, Ginga 
and RXTE data (Apparao et al. 1991; Yamauchi et al. 
1990; Hall et al. 2000). However, when compared to those 
previous observations, we find a larger value for Ah- This 
difference can be explained by the fact that the A^h mea- 
surement obtained with BeppoSAX is substantially more 
accurate thanks to the better sensitivity and spectral cov- 
erage at low energies afforded by the BeppoSAX LECS 
and MECS. Thus, in our opinion the A^h value obtained 
with BeppoSAX should be considered as the correct one. 

We also note that the Nn values we found in both 
flare and low-level spectra are substantially higher than 
the Galactic one along the 4U 0114+65 line of sight, which 
is Nf[ — 0.76xl0 22 cm~ 2 according to Dickey & Lockman 
(1990). This fact is also confirmed by the optical V-band 
absorption estimate along the LSI +65°010 line of sight 
made by Reig et al. (1996): the value they obtain, Ay = 
3.84, corresponds to a Afa = 0.69xl0 22 cm~ 2 using the 
empirical formula by Predehl & Schmitt (1995). This is 
consistent with the Galactic Nh value and substantially 
smaller than the Ah we found from X-ray spectral fits. 
This difference is very likely due to the fact that the stellar 
wind absorbs soft X-rays, but does not cause an increase 
in E(B — V), and thus Ay, because the wind does not 
contain dust (see, e.g., Cox et al. 2005). This fact further 
suggests the presence of a noticeable and variable amount 
of neutral gas local to the system and very close to the NS 
(possibly, the accretion stream itself). 
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A variation in the characteristic energy E cn t and (to 
a lesser extent) in the values of -Efoid an d T of the 
HighECutPL model is also found, with these parame- 
ters decreasing in value as the total flux decreases. This 
trend is opposite to what was found by Hall et al. (2000) 
in the RXTE data. This may at least partly be due to 
the wider spectral coverage afforded by BeppoSAX, which 
allowed a better sampling of the X-ray spectral states of 
4U 0114+65. 

Concerning the possible presence of a further compo- 
nent in the softer part of the spectrum (below 3 keV), 
we note that Hickox et al. (2004) found that a soft ex- 
cess, usually modeled using a BB with feTee ~ 0.1 keV, 
is ubiquitous in X-ray pulsars (XRPs) and, in particu- 
lar, in HMXBs hosting pulsating NSs. According to them, 
XRPs with luminosities larger than ~T0 38 erg s _1 are 
produced by reprocessing of hard X-rays emitted by the 
NS by optically thick accreting material (e.g., an accre- 
tion disk) , whereas in cases of XRPs with luminosities be- 
low ~10 36 erg s _1 this mechanism can be excluded and, 
most likely, the soft excess arises from different processes, 
namely emission from either a photoionizcd, or collision- 
ally heated, diffuse gas or from the NS surface. 

Given the size, exceeding the typical NS radius of a 
factor ~20, that we found for the thermal component re- 
sponsible for the possible soft excess in 4U 0114+65, we 
suggest that a cloud of diffuse plasma around the NS pro- 
duces this excess. 

The Fe emission, its variability and the fact that it 
is better seen during low emission levels were noted by 
Yamauchi et al. (1990) and Hall et al. (2000). This emis- 
sion line may arise from the heated gas cloud depicted 
above as the possible cause of the soft X-ray excess and, 
as with the continuum emission from this cloud, is ap- 
parently better detected when 4U 0114+65 is in the low 
state. 

We did not find a CRF at 22 keV in any of the two 
brightness states; this is at variance with the findings of 
Bonning & Falanga (2005), who do not report any value 
for the line depth of their CRF detection. Therefore we 
are not able to confirm the estimate of the NS magnetic 
field strength made by these authors and, in turn, the 
magnetar origin proposed by Li & van den Heuvel (1999) 
for the NS hosted in this system. 

From visual inspection of the right panel of Fig. 4, 
there is a suggestion of a feature at ~20 keV in the low- 
level spectrum of 4U 0114+65. Using the NCR technique, 
we demonstrated in Sect. 3.3 that this feature is not due 
to an incorrect modeling of the X-ray spectral continuum, 
and thus we exclude the presence of any absorption feature 
at ^22 keV in our BeppoSAX data at both flux levels. 

This same situation is observed for another X-ray bi- 
nary pulsar: Vela X-l. As pointed out by Orlandini & Dal 
Fiume (2001) and Coburn et al. (2002), in the 15-30 keV 
energy interval there seems to occur a change of slope that 
is not well described in terms of actual spectral models, 
and that can give rise to residuals that can be misinter- 
preted as absorption features (see also the discussion on 



the X-ray binary pulsar OAO 1657—415 in Orlandini et 
al. 1999). 

The above findings suggest the following scenario for 
4U 0114+65: the accretion stream from the secondary 
star, in the form of a low angular momentum stellar wind 
as suggested by, e.g., Reig et al. (1996), is captured by the 
accreting NS and creates a plasma cloud around the ac- 
cretor. This cloud thermally emits at a temperature kT ~ 
0.3 keV and also produces an iron emission line. The ther- 
mal emission of the cloud is best seen at low hard X-ray 
flux levels, when this plasma cloud is denser, otherwise 
it is overwhelmed by the more intense underlying harder 
X-ray emission from the NS surface and/or from the NS- 
cloud boundary layer. 

The total luminosity of the source is then driven by 
the density of the accreting cloud: when higher, it has the 
double effect of masking the inner harder emission and 
of producing a more intense softer thermal emission from 
its outer parts. Within this scenario, the presence of the 
superorbital period reported in X-rays by Farrell et al. 
(2005) cannot be explained by the presence of a precessing 
warped accretion disk; rather, it might be produced by the 
NS spin axis precession. 
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